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The heat transfer characteristics of a rocket motor 
storage container system have been investigated using 
Eu vtIcal and experimental techniques.  Analytically, both 
closed form and numerical solutions have been developed. 
These solutions may be used to determine maximum tempera- 
wres and temperature gradients within the rocket motor. 
Comparison between theoretical and experimental values of 
temperature are within the estimated experimental uncer- 
ο Ετος ος t 3° F. It is proposed that the theoretical 
Solutions can be used to thermally optimize container 
design. 

A secondary investigation was carried out to determine 
Ame easib1lity ok uslag cholesteric liquid crystals, a 
Amper cure sensitive material, to thermally map the surface 
of the container. The crystals appear to remain stable 
under desert type conditions and produce brilliantly colored 


dis pillas of the temperature field., 
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eee RODUCT ION 


The purpose of this investigation was to develop a heat 
transfer model that will allow prediction of the temperature 
Mestribution in a container stored rocket motor placed in a 
hostile thermal environment such as the desert. It is pro- 
posed that such a model would be a useful tool for thermally 
Setrmizing future container designs. As extreme variations 
in the rocket motor temperature may lead to large thermal 
Stresses in the propellant which could result in fracture, 
otherwise degrade the performance of the motor, a major 
objective of this study vas to design a model that could 
mora billy predict the thermal gradient in the motor. The 
predictions would be based on the surface temperature dis- 
put on, the thermal properties and the geometrical de- 
tails of the system. The model may also be used to predict 
a critical temperature range over which the propellant must 
be chemically stable when in a storage situation. The 
perimit Of thus temperature range is referred to as the 
design temperature of the system. As the design temperature 
for most weapon development projects is derived from dump 
storage COR SO AS. a dump storage situatíon was used to ob- 
ΕΤΕ experimental data for this project. 

Several approaches were taken to predict the rocket 
πι ο τ stcemperaturesadresrsbution from a knowledge of only the 
surface temperature distribution of the storage container 


andsche thermal properties and geometrical details of the 


TI 





experimental model. The experimental model used in this 

Nest was a wird Navy antisubmarine rocket (ASROC) 
motor, filled with dry desert blow sand to simulate the pro- 
Alant, and placed in its storage container. This container 
System was placed in a dump storage site at the Naval Weapons 
Memrer, China Lake, California to simulate a desert environ- 
went. 

The method of complex temperatures [Ref. 1] was used to 
Gevelop an analytical prediction of the transient tempera- 
οι Ἱσοσιὰά that exists in a container stored rocket motor. 
Weemanalytical model assumes that heat is transferred only 
mer radial direction and that the container surface tem- 
Perature variation is sinusoidal with time. Comparison 
between theory and experiment is within experimental uncer- 
Bu when temperature is interpreted as "bulk" temperature. 
miemanalytical model is especially useful for studying 
ne .5:1 and thermal physical property effects on rocket 
motor temperature. Such parameter studies have been carried 
ande resultes are presented in a form that will be 
useful from a container design point of view. 

ΠΤ, :, wand ot, a computer program for transient 
and steady-state temperature distributions in multidimensional 
ο πο σου used to obtain detailed information about the 
ο πο εοες ο the rocket motor. TRUMP allows actual con- 
tainer surface temperature distributions to be used as well 
as sinusoidal variations. In addition, both one dimensional 


cuu I9udimenstonal (radial and circumferential) 


T2 





Mt ranster were modeled with TRUMP, using both the sinu- 
wea t and actual temperature distributions. The actual 

Aa erature distributions were obtained from the experimental 
Aia of the motor EU en system. 

Comparisons between the experimental values and those 
predicted by the models were in good agreement ,with those 
τες by TRUMP using the actual temperature distribution 
as the boundary condition being the closest. However, the 
EMENUSOvdal variations used in both the analytical MEET and 
BNSSUPRUMP model are also suitable for design purposes. 

Another aspect of this project was to obtain the storage 
emtazner surface temperature distribution using cholesteric 
τη crystals, a material that undergoes brilliant changes 
ICO lor over known, well defined temperature ranges. Color 
Slides and movies were taken of the liquid crystals demon- 
ciao the feasibility of using them for on site tempera- 


ture measurements. 


13 





: 


LL. BACKGROUND 


Πο ο the Naval Weapons Center, China Lake recognized 
the need for a concerted attack on the problem of thermal 
Eme la assienment for new weapon systems. In 1963 a task 
force was established to study the complete environmental 
ΙΤ. ia determination problem. The key to this problem 
Seema to be the thermal area in the storage and transporta- 
events of any item. It was realized that transportation 
was a short term situation compared to the storage situation. 
Therefore, the major portion of the life of an item must be 
storage. There are three types of storage; covered, igloo 
and dump. The dump storage situation leads to the more 
extreme thermal exposure situations which then leads to the 
design temperature. 

AS data was not available for the dump storage situation, 
instrumented storage dumps were created at representative 
pu conva worldwide basis so that statistical data could 
Ῥε ο πνεά on a variety of ordnance. The first site was at 
China Πίο. California, in the middle of the Mojave Desert. 
This site now has the capability to return about 250 channels 
of information on a continuous time-temperature basis (Fig- 
it Other arctic and tropical sites were set up to study 
exbemaedae Conditions: 

The dump storage situation was reproduced to study the 
extreme situation. The ordnance was exposed singly, directly 


situated on the ground, with the long axis aligned in the 
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north-south direction to allow maximum normal exposure of the 
Ata ner surface to the sun s rays. In actual practice, 
Smemance zs usually stacked and oriented in other than a 
north-south ο ο e Tebo ao ene extremne situa- 
tion. Ordnance sitting on the ground receives reflected 
Bdrarron from the ground, cannot quickly give off heat by 
Eouuductrion to the soil, and is not as apt to be cooled by 

pre prevailing breeze; therefore, extreme temperatures 
result. 

The most important source of heat to the ordnance is the 
Ber radiation from the sun, with reflected radiation of 
ο απσαιτν importance. For extreme conditions to occur the 
να πππε be calm (less than 5 knots), the sky clear, and 
the outside air temperature high. After sunrise, the ord- 
nance skin temperature rises much more rapidly than the 
MOC nt air temperature; therefore, the surrounding air 
eeollae =the ordnance, rather than heats it. 

The rocket motors used for the tests were military sur- 
peus Even though the material had served its intended in- 
Fleet purpose, it was still representative of new hardware, 
when viewed in a thermodynamic context. When inert rocket 
motors were available, they were used intact; however, in 
most cases, once-fired hardware was used. Thoroughly dried 
desert blown sand, being similar in thermal properties to 
most propellants, was used to backfill empty rocket motors. 
It was assumed that the thermal response of the sand filled 
Motors was essentially the same as actual propellant filled 


motors. 
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Piece MENTAL PROCEDURE 


Although Naval Weapons Center, China Lake had accumulated 
Vast amounts of data in the past, it was decided to instru- 
ment a rocket motor storage container system especially cn 
this project. This would allow base data to be taken exactly 
where it was required. It also allowed variations in the 
System without interfering with one of China Lake's ongoing 
eects. An ASROC system was chosen for this study. The 
outer storage container was 75 inches long vith an inner 
gm erter of 18 inches and a wall thickness of 1/16 inch. The 
rocket motor was 57 inches Jong with an outside diameter of 
ches and a wall thickness of 1/4 inch. Both the con- 
tainer and motor were made of steel, 

The rocket motor storage container system was instrumented 
with 20 gage copper-constantan insulated thermocouple wire 
πας απ ISA Calibration of + 1-1/2°F over the range -75 
to +200°F. Twenty-one thermocouples were originally placed 
emsther system with positions indicated in Figures 2 and 3. 
The ambient air temperature was meaSured with thermocouple 
number 19 which was located in a Stevenson shelter about 60 
feet a av from the system (Figure 4). The thermocouples were 
mounted intrinsically on the motor and storage container. 
-Two small holes were drilled approximately 1/8 inch apart in 
the metal and the individual wires were inserted in the holes. 
The metal was then hammered around the wires until a snug 


fit was obtained. Bead thermocouples were mounted at the 


dy 





κ L 
2 0° pe 
2 
Steel > e 
7 
INO 22 1139 30 9 ES 
Sand 
E 
9 
Air Gap 
10 
Figure 2. Thermocouple Locations on Experimental System. 
Five thermocouples were located under the section painted 
with the liquid crystals. Their locations corresponding to 
the ones shown above are: #14= #1, #15= #2, #16= #8, #17= #9, 
and #18= #3 (See Figure 3). 
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Member of the motor and in the air gap. The thermocouples 
located at the center of the motor were supported by small 
Peces of wood several inches from the head. The use of 
these supports was necessary to keep the thermocouples in 
position when the motor was being filled with sand. After 
all the thermocouples on the rocket motor were in place, 
ο et MoLOnewas filled with dry desert blown sand. 

The wires from the two thermocouples located in the center 
15 motor were led out a hole in the end cap. To avoid 
Settling of the sand after the motor vas in place on the 
Sica a resulting air gap being formed between the 
sand and the motor skin, the sand was compacted by striking 
the sides of the motor with small sledge hammers and then 
ae additional sand through the hole in the end cap. 
This was continued until the sand was tightly packed, The 
ντι the end cap was then sealed. The rocket motor was 
Comemully placed in its storage container (Figure 5) which 
had previcusly been instrumented with thermocouples. The 
thermocouples in the air gap were mounted a affixing the 
lead wire to the rocket motor at the desired position and 
then putting a 90 degree bend in the wire so that it placed 
the bead of the thermocouple approximately 1.5 inches into 
the air gap. Neither the thermocouples in the center of 
the motor nor those in the air gap could be considered accu- 
rately positioned; however, every effort was made to mini- 
mize positioning errors. All thermocouple wires were 


located inside the storage container and were led through a 
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Maletin one end. This hole was then sealed. The two halves 
Seethe storage container were then bolted shut. 

Mmemeouter surface of the rocket motor and the inner and 
outer surfaces of the storage container were all painted 
various shades of haze gray. Weathering had caused the 
wonted surfaces to appear fairly rough. This is typical 
urhe conditions of a storage dump. From the condition of 
the surfaces, it was estimated that the emissivity was 
approximately 0.9. 

Prion, to loading the rocket motor into the storage con- 
tainer, it was decided to apply liquid crystals (See Appen- 
WRA to part of the storage container surface in order to 
obtain a thermal mapping of the surface temperature at any 
instant of time, Liquid crystals are temperature sensitive 
materials that produce immediate thermal images in a pattern 
of colors which respond rapidly to minute changes in sub- 
strate surface temperatures. A second reason for applying 
Meery stals to the container surface was to determine the 
feasibility of using the crystals under adverse environmen- 
d onditions (desert atmosphere). Prior to applying the 
als ao ap ch strip of the storage container, 20 inches 
from one end, was sprayed with two coats of Testors Spray 
Pla Enamel No. 1249, Flat Black as a background for the 
2] one inch strip of ll different ranges of crys- 
τ... .15; αρρτοπίπαεειῖν 1/2 inch of black paint between 
men was applied over the black paint. Two coats of each 
crystal were applied, using a small paint brush. The first 


coimas allowed to dry completely before the secona coat 
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LE asprplied. After the crystals were dry, two coats of Rez 
Bemyurethane (gloss clear plastic coating, interior-exterior 
77-5) coating were applied by brush completely covering the 
crystals and black painted area. The polyurethane coating 
weomapplied to protect the crystals from wind blown sand and 
from the ultraviolet rays of the sun. Ten of the eleven 
Eis had been previously calibrated [Ref. 41]. Using the 
Comstant temperature bath procedure recommended in Ref. 4, 
as Calibrated and the complete calibration results are 
shown in Table I. 

The rocket motor storage container system was then moved 
to the China Lake dump storage site. The system was aligned 
Ino rth=south direction, well away from the influence of 
Sewer ordnance (Figure 6). The thermocouplée leads were 
connected through a junction box and underground cable to a 
ome” ell Electronik 25 Recorder which had been calibrated 
to read the thermocouple output directly in degrees Fahren- 
Dto an aecuracy otf 4^ l'"F. The recorder was located in 
an air-conditioned building about 60 feet from the system. 

Initial data indicated that the number 7 thermocouple 
Was not responding properly and therefore this data was 
P ceda Initial color photographs were. taken of the 
liquid crystals and it was immediately apparent that good 
dma mappings could be obtained if the crystals were 
stable under the adverse desert environment. The brilliance 
of the colors exhibited by the crystals under the bright 


desert sun was much better than had been expected. The 
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τον αι 
Calibration Of Liquid Crystals 
Calıbration Bath 


NCR Manulacturer s ye Goats 
Desig. Color Change Responses quita Crystals 
Ec RG 
R-27 | 
Red 200 2009.5 
Green Da ο 28.0#.5 
Blue 50.0 28.7+.5 
R-33 
Red 393.0 S247 44) 
Green 34.6 33.3+.5 
Blue 36.0 34. 24-. 5 
Β-37 
Red 31:50 36. 2. 5 
Green 38.6 1.» 
Blue 40,0 38.0+,5 
R-41 
Red απ ο ο -ᾱ 5 
Green 42.6 41.0+.5 
δη ιο 44.0 AOS 
R-45 
Red 45.0 42.8+.5 
Green ΠΟ Ὁ OS 
Blue 48.0 44.3+.5 
R-49 
Red 49.0 ΠΕ 71.5 
Green 50.6 47.1+.5 
Blue 52.0 48.4+.5 
R-53 
Red 5370 mes ᾗ 
Green δ 5 S201» 
Blue 5670 ΠΠ. 
R-56 
Red 56.0 53.8+.5 
Creen 5:25:06 3620.25 
Blue 59.0 DO Got. 
R-59 
Red 59.0 4.5 
Green 60.6 ISF 
Blue 62.0 Oto i 5 
S-62 
Red 62.0 ου τη. 5 
Green 6 60.4.5 
Blue 63.0 SE o5 
S-64 
Red 64.0 60.9+.5 
Green 64.6 61.4+.5 
Blue 65.0 6ος 
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Experimental System at Dump Storage Site. 


Figure 6. 





system was allowed two weeks to reach a periodic steady 
ο τε before additional photographic data was obtained. 

Extensive photographic data was collected on 27 and 28 
July 19/72 after two weeks of exposure to the desert environ- 
E... Both super 8 mm and 16 mm color movies and 35 mm color 
peemdes were taken of the liquid crystals. No colored filters 
were used on any of the cameras, although standard haze fil- 
ters were used to take the super 8 mm movies and most of 
the 35 mm slides. 

AS times a second storage container, ee with- 
οσπι αᾱ rocket motor inside, was instrumented with intrinsic 
Ehermocouples in the same manner as the previous container. 
As uv three data channels remained open on the recorder, 
only three thermocouples were applied to this new container. 
The three thermocouples were applied at the 0300, 0900, and 
1200 positions at the midpoint of the container. ntsc on 
tainer was set end to end with the system that was already 
in place at the site. The purpose of this study was to 
Humcrmrirne if the inclusion of the rocket motor in the con- 

τ: παπα a significant effect on the surface temperature 

of the container. Thermocouple #7 was connected at the 0900 
mendon, 723 at the 1200 position, and #24 at the 0300 posi- 
tion. It was immediately apparent that thermocouple #7 was 
emning to give unreliable readings and therefore the 


data taken on channel #7 was neglected. 
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ο ου ΟΕ ΤΙΝΑ ANALYSIS 


A. ONE-DIMENSIONAL ANALYTICAL MODEL 

The first step was to try to devise an analytical model 
that ια simulate the actual rocket motor storage contain- 
er experimental system. The first simplifying assumption 
was that the storage container temperature could be modeled 
by a sine wave which had a period of 24 hours. A compari- 
nof the sinusoidal variation to the average (bulk) stor- 
age euer temperature [obtained by averaging the four 
Srermocouple readings on the surface of the container (1, 
LOS and 12) as shown in Appendix D] is given in Figure 7. 

The method of complex temperature as presented by Arpaci 
EC vas used to find the steady periodic solution of a 
body experiencing a periodic sinusoidal disturbance. A com- 
plete analytical derivation is given in Appendix B. The 
general heat conduction equation in cylindrical coordinates 
ΤΠΕ basis for this derivation. It was assumed that 
there was one dimensional radial heat flow with no conduction 
Meee axial or circumferential directions, that no heat 
sources existed in the model, that the rocket motor storage 
Container system was infinitely long, and that the sinusoid- 
ally varying surface temperature was spacially uniform over 
the entire container surface. The storage container tempera- 
ture is assumed to vary as 


- E i + 
DS as a) sin wt T 


A 


28 





-aıny3easadusa] ATNY 02 UOTIBTARA SinjB19duSÍ TEPprosnufs 340 uosraeduoo 


(SAINUTK) SUTIL 


007 T 0OZI 0007 008 009 007 ΠΣ 
ο 
ο 
ο το ο 
ο 
ου 
s os 
ο 
ο No o, 
ο 
ο 
ο 
e ο 
ο 

a ο 

: ο 

5 ο 

Ὸ 
ο 
ο 
2 ο 
ο ο 9 
ZLZ6T 3Sn3ny Z FO 5 9 
o” e ο 


eze( sin3e19aduaz] yng 5 B 


oinjgezeduep [9?prosnurg——— 


“/ 9anarı 


06 


OTI 


OET 


(1,) eanaezedue] 


29 





where Τμ = maximum bulk temperature of the storage 
Container 


Ty ο ο ο οσα temperature of the storage 
container 

Oee=strequency of the sinusoidal variation 
(21/24 hours) 

t = time (hours) 


It was assumed that all the thermal properties remained 
constant over the temperature range of the problem. The 
effective heat transfer coefficient, h, across the air gap 
between the storage container and the rocket motor combines 
Mhe heat transfer effects of radiation, convection, and con- 
1 into one coefficient. The radiation coefficient was 


linearized by assuming constant representative temperatures 


in the equation 


2 2 
= E + 
PRAD 39 DE απ M 
where 0 is the Stefan-Boltzmann constant and Ps LS ele 
radiation êxchange factor. The convection coefficient is 


Ene effective conductivity if ee obtained from the eee 
mannteorrelations [Ref. 5], divided by the width of the air 
Pape in the analytical model, the effective conductivity 
was assumed to equal the conductivity, thereby treating it 
as pure conduction and giving the equation 


ON 


Aca condition was not specified in this deriva- 


enas the only concern was the steady-state, periodic be- 


havior. The steady-state solution is (Appendix B) 


T(r,t)-T, YBER * (a£) *BEi * (a&) 
Ə (r,t) = ETT ER c sin(wt+6*) 
== A 
M A Xp + X, 
Meee) =O sin (at + 6*) 
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where T(r,t) enea perature of a point r in the rocket 


MOL Or Tat time t 





wr 3 
a = d = conduction parameter 
É uM 
E = ie dimensionless distance from the center 
o of the rocket motor 
a aoe radiuso the rocket motor 
e el ance from the center of the rocket motor 


α oi tbernalgsd3sffusivity 
p = density 

k = thermal conductivity 

e 


= specific heat 











BER = real Bessel Function 
BEi = imaginary Bessel Function 
a a 
== 4 4- : . 
Χρ BER (a) BER, (a) + E? BEi, (a) 
X. = BEi (a) + A Rea) 
i O 1 2p 1 
hr 
B = τ = Biot modulus 
-1 BEi (a£)X, ~- BER (a8) xX, 
Grete 


BER (26) Xp F ΡΕ1 (αξ)Χ, 


Two computer studies were done based on the steady state 
solution. The first study was a completely dimensionless 
Situation which served aS a parameter study of the effects 
Gmevaryine a and B on the temperature and the time lag of 


the temperature at various positions in the model. 


- 








a = ^ua = conduction parameter 
and n 
hr, 
B = E = Biot modulus 
Parameter a was varied from 1.0 to 5.0 and ß was varied 
from 0.1 to 100. These were the only values studied, as 


BI 





only values within this range are of interest in this type 
Problem. The computer program and its output are given at 
ο επά of Appendix B. The output lists the following 
Values: 

1) a, the conduction parameter 

2) RB, the Biot modulus 


3) &, the non-dimensional distance from the center of 
the motor 


4) 6%,the time delay between the maximum storage con- 
parmeratemperature and the maximum temperature 
ΓΓΕΤ“ point of interest in the motor 

5) 9 ,the relative amplitude of the maximum tempera- 
ture at the point of interest compared to the 
maximum temperature of the storage container 

The time delay is given in radians, where 2T radians equals 
one complete cycle. A graph of the time delay versus P for 
ENconstant value of "a" is given in Figure 8 at three dif- 
trent positions within the motor. A graph of the relative 
amplitudes of the temperatures versus ß for a constant value 
as given im Figure 9. It was noted that for a con- 
stant value of "a", the time delay decreased as B became 
larger. As the point of interest approaches the center of 
the rocket motor, the time delay increases. The relative 
amplitude of the temperatures also becomes larger as ß is 
mMercased when the value of "a" is held constant. If Bis 
held constant and "a" is varied, the time delay increases 
and the relative amplitude decreases as "a" increases. 


The second study was obtaining the analytical solution 


to the particular rocket motor storage container system 
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died at China Lake. The thermodynamic properties of dry 


End were obtained from Ref. 6 as 


MN E NIE. 
k = 0.188 BTU/hr ft °F 
c = 0.195 BTU/lbn °F 


Substituting these values and using 1440 minutes (24 hours) 


as a complete cycle, the parameters a and ß were calcu- 











lated for this model as 
Or, 
a = 22003 
a 
where τσ 5.75 inches, the inner radius of the rocket motor. 
T 
B = PR σι 
where h = h coN a BRAD 
k 
and hogy * q * 6.48 x οσα — 
: hr-ft^*?F 
uere fr = 2.94 inches, the distance across the air gap 
-2 BTU 
and KAIR = 1,62 x 10 OTE 
> 2 2  BTU 
ο otros TUR CET, s 1.09 B 
In r^ FE 


where O is the Stefan-Boltzmann constant, UN is the 


radiation exchange factor which for this geometry is 


= 0.84 
- 1) 


when Ec &> = =: re 60, rj - DIU 


therefore h - 1.15 BTU/hr ft??F 
The average surface temperature of the storage container 


No un dto 5 N04 F for a particular day at China Lake, 
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with a maximum temperature of τ» F. These values were ob- 
Med by averaging the readings of thermocouples 1, 8, 10, 
and 12 as shown in Appendix D which give the bulk tempera- 
Eure. | 

mime temperatures of seven positions within the rocket 
Motor were calculated and the results are printed at 30 
mute intervals for one complete cycle in Appendix B. A 
graph of temperature versus time was plotted by the computer 
Bnowine the relationship between the surface temperature of 
the storage container (TINF), the temperature on the outer 
ΜΠΕ the rocket motor (TEDG), and the temperature at the 


meme of the motor (TCEN). This graph is Figure 10. 


Pee TRUMP MODEL 

Buccrocket motor storage container system at China Lake 
was modeled ou TRUMP, a numerical conduction code, (See 
Erondrx C for a description of the TRUMP program) to pre- 
dict the temperature at any point in the system from a know- 
ledge of the storage container surface temperature variation, 
Enesthermal properties and the geometrical details of the 
system. Two models were used to simulate the rocket motor 
storage container System and several variations of each 
model were investigated. 

The first model assumed one dimensional heat transfer 
radial). The system was modeled as two infinitely long 
Concentric Steel cylinders, the inner of which was filled 
e dry sand, A 2.94 inch air gap separated the cylinders. 


The model was subdivided into concentric volumetric elements 
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with representative nodal points as given in Figure 27, 
appendix C. It was assumed that the storage container sur- 
meee temperature was spacially uniform. From the data given 
in Appendix D and the observation of the liquid crystals' 
thermal mapping, it was obvious that the temperature distri- 
Eon on the storage container vas not spacially uniform. 
order to »simulate a spacially uniform condition, the 
readings of the thermocouples located at the 1200, 0300, 
wood, and 0900 positions (#1, 8, 10, and 12) were averaged 
and this average value of the surface temperature (referred 
to as the bulk temperature) was used as the spacially uniform 
Benperature distribution. Two methods were used to describe 
the container temperature. The first method used the maxi- 
mum bulk temperature (138°F) and the average bulk tempera- 
Pere (104°F) of the storage container to generate a sine 
Wave with a period of 24 hours (1440 minutes). The second 
metned took the bulk temperature readings at two hour inter- 
dals and fed this data into the TRUMP program in a tabular 
(temperature versus time) form. The version of TRUMP used 
in this problem was limited to a table length of 12 tabular 
values. TRUMP interpolated between the tabular points. 
Figure ti compares the actual bulk data with the sinusoidal 
approximation and the interpolated tabular values. 

several assumptions were made to simplify the solution 
Athis problem. As the thermocouple data from the storage 
container gave an average value of the temperature across 


the 1/16 inch steel wall, node 12 was modeled as a zero 


37 





*1030U 2143 
JO 193Uu92 əy3 39 einjq?ieduoj eua ST NuolL PUE Cc9030099939029919 570 920135 319d 99 
s9eyrıns 943 ST DUAL “19UTEJUOD eo2e103S 5913 JO 31n3et19du93 soejins 8913 ST ANIL 

atoum *9urry u3r4 VUOTIPTABA DANJEISAUDL JO UOTITPOIA TEITSZATEUY TL Ot οαυσοτα 


ο Τ DOSE 0001 | 008 009 007 007 


NIJI 2 


Ὀ 31 


ANIL 


06 OL 
(ας) eanjerodusy, 


OTT 


OET 


38 





* SUOTIRUTXOAAdY ¿INNÁL OMIL oa oina?aoduoep xIng JO UOSTARPdUO) 


(SIINUTHN) SUTIL 


00*I QUSE 0001 008 009 00% 
ο 
ο 
ο To 
o Y ο. 
Ζ o 
c, “S 
/ ao g 
No 
* ο 
Ao 
/ \o 
Y Χο 
\ 

\o 
\ 
Q 
\ 

\ 

ολ 

ο 
\ 
9 
\ 
uotjzeurxoiddy ıeTngeL ---- N 
\ 
oo 
uorzeutTxoiddy [eprosnurg -——— SAC 


e}eq [?3uosuradodxg o 


ST Sp Ed 


00 


C 


DUETO 06 08 04 
(1,) »ınyeıoduaL 


OYT 


39 





volume boundary node with a known temperature impressed upon 


mee Lt was also assumed that heat transfer across the air 
gap occurred by radiation and conduction alone. Free convec- 
tion effects were initially neglected. This assumption was 


Ie modified to investigate the free convection effects. 
All surfaces of the storage container and the outside surface 
πο rocket motor were painted various shades of haze gray 
and it was estimated that the emissivity of these surfaces 


was 0.9. The radiation exchange factor, | por chis 


p 
model was the same as that for the analytical solution E 
= 0.84). It was also assumed that there was perfect thermal 
contact between the rocket motor and the sand that filled 

w This neglects the possibility that the sand might 
ENntiy settle after being on the site for a long period of 
time. 

The second model assumed two dimensional heat transfer 
51 and circumferential). The same physical model was 
used as in the one dimensional case with the sole exception 
that 48 nodes were used instead of 12. The representative 
nodal points and an example of the thermal connections from 
one of. the nodal points are shown in Figure 28 in Appendix C. 
The four nodes on the surface of the storage πον νοχὰ 
modeled as zero volume boundary nodes. The sinusoidal and 
tabular representations were used to describe the surface 
temperature of the storage container at each boundary node. 
wnal data taken at each position, rather than bulk data, 


were used as the input data for these representations. 
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The same assumptions made in the one dimensional case 
were also applicable to the two dimensional model. A com- 
plete discussion of the calculation of the radiation ex- 
puunmpe factors in the Ho dimensional case is given in 
Appendix C. 

The effect of natural convection was studied in both 
the one end two dimensional models. References 5 and 7 
give correlations between the Grashof number based on the 
gap width and the effective thermal conductivity. The 
Grashof number was calculated from the equation 


p* gB(AT) 6° 


Gr = u? 


where ὃ Wee home tive air gap 


AT = maximum temperature difference at any instant 
of time in the air gap 
m o 
το. where T 9 565 R 


T 


EET — 565"R, air has the following properties 


0.07 mes 


M 


p 


u 0.046 lbm/hr-ft 


The maximum Grashof number for this experiment was calcula- 


6 


ο το be 1,25 x 10 . The diameter ratio was approximately 
Sand the log Gr = 6.1. From the Beckmann correlation 
beet. 5]. the effective thermal conductivity ratio ES) was 
Ao imately 3,2, Using the Liu correlation [Ref. 7] 

ae er κ y E 

k ` l.36+Pr i 


where the Prandtl Number = 0.707. An effective thermal con- 


ductivity of 4.0 was assumed as the average value of these 
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myo, correlations and it was used to study the effects of 
free convection. This change was placed into the TRUMP pro- 
meame by increasing the value of the thermal conductivity of 


EN Nby a factor of 4 in each of the TRUMP runs. 
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V. d RESULTS 


fee ANALYTICAL MODEL 

Using the sinusoidal temperature distribution as an ap- 
proximation to the actual average experimental data as 
Shown in Figure 7, comparisons were made between predicted 
temperatures and actual temperatures for two radial loca- 
Bons in the rocket motor. Figure 12 compares the results 
ΙΕ surface of the rocket motor and Figure 13 does the 
comparison at the center of the rocket motor. An uncertain- 
ty analysis is given in Appendix E which establishes the 
Epucertainty bounds for both the predicted and the actual 
temperatures. These uncertainty bounds are included in 
Figures 12 and 13. 

It is readily seen from Figure 7 that a sine wave was 
not an ideal fit as an approximation to the experimental 
[mane as it varies as much as 20°F during part of the cycle. 
However, it was also noted that the sine wave closely en 
ped the experimental data during the heating phase of the 
cycle and only during the cooling phase were there large 
variations. As the main purpose of this study was to design 
a model that would be useful in optimizing storage container 
design, the errors in the cooling phase are not critical as 
long as the temperatures reach the same minimum point before 
beginning another cycle. Figure 12 shows that the maximum 
Burrace temperature of the rocket notor predicted by the 


analytical model is a good approximation to the actual 
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E rental data. Again it is noted that, in actuality, 
the motor cools faster than the predicted value. The maxi- 
E difference in temperatures on the surface of the motor 
ms 15°F. Figure 13 shows that the predicted value and the 
experimental value of the temperature at the center of the 
Bucket motor were in close approximation except during the 
early stages of the cooling phase where a maximum tempera- 
Nc variation of about 5°F occurred. 

One of the reasons the system cools faster than predicted 
could be the light breeze that is usually evident in the 
early afternoon hours at China Lake that is not present dur- 
ing the morning. No attempt was made to shield the system 
from the wind to study the effects of a light breeze on the 
Surface temperature of the storage container. 

Ao ter point not taken into account by the analytical 
Model is the fact that the time deley at any point in the 
system is not constant throughout the day as predicted in 
Figure 10, but varies as given by the data in Appendix Di 
Time delays between the peak temperature on the container 
Surface and the peak temperature at the center of the rocket 
Motor vary from about 250 to 400 minutes, whereas the lov 
temperature on the surface of the container and the low ten- 
Ar ature at the center of the rocket motor vary from about 
HO to 250 Minutes. The analytical model predicts a constant 
An iraton Of 360 minutes at the center of the rocket motor 


ema 159 minutes at the surface. 
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B. TRUMP MODEL 
1. One Dimensional 

Four variations of the one dimensional TRUMP model 
were investigated and compared to the experimental data. 
Figures 14 and 15 compare the TRUMP predictions to the act- 
meee xperimental data at the surface and the center of the 
Pocket motor,respectively. The TRUMP variation used for 
this comparison modeled the storage container temperature 
with tabular data (See Figure 11) and assumed convection 
was present EE. foe The uncertainty analysis (Appendix 
established the uncertainty bounds for both the experi- 
eemal and the analytical data in these Figures. The varia- 
EMnSbetween the bulk temperature predicted by TRUMP and 
the experimental data closely matches with only two experi- 
meade points in Figure 14 falling outside the uncertainty 
bounds for this one dimensional model. Figure 11 shows 
meee the tabular data that TRUMP interpolates is a good 
EDproximation to the averaged experimental data. At the. 
center of the motor, as shown in Figure 15, all experimental 
Points fall within the predicted error bounds. A comparison 
ethe four one-dimensional TRUMP variations are given in 
Ex cures J6ócand 17 at the surface and the center of the rocket 
Motor respectively. It is clearly seen from these Figures 
Ehat the convection assumption results in an increase of 2°F 
in the Maximum temperature and a decrease of 2°F in the 
minimum temperature on the surface of the rocket motor. This 


Bonperaturescnange drops to # 1.5°F at the center of the 
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motor as shown in Figure 1/. The differences between the 
sinusoidal approximation and the tabular approximation of 
the experimental data was clearly shown in Figure 11, The 
mired in Figures l6 and 17 can be easily correlated to that 
in Figure 11, thereby explaining the differences in the pre- 
eted values. 
2, Two Dimenstomel 

Eonnmsevarıatıons of the two dimensional TRUMP model 
were investigated and compared to the experimental data. 
Momparisons of each TRUMP variation to the experimental data 
are given in Figures 18 and 19 for node 8 (located on the 
τη of the rocket motor at the 1200 position) and node i 
(at the center of the rocket motor) respectively. These 
Figures show that the TRUMP variations that used tabular 
data to model the surface temperature of the storage con- 
ἵατπες predicted temperatures that more closely approximated 
the experimental values than were those predicted by TRUMP 
Mar tations using sinusoidal data to model the surface tem- 
perature. Appendix D shows that all points on the surface 
of the storage container reach their minimum temperature at 
the same time; however, these points reach their maximum 
temperature as much as 200 minutes apart. Whereas, ali the 
points on the surface of the storage container are in phase 
at the minimum temperature, they rapidly become out of phase 
as the container temperature rises. This varying phase 
Shif makes it difficult to model the four boundary nodes 


with sinusoidal approximations which must have constant 
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phase shifts. Sizable errors in the input data during some 
parts of the cycle were caused by these varying phase shifts. 
Messe errors in the input data led to the variations in the 
predicted temperature values. As noted in the one dimen- 
Sional section, the inclusion of convection effects does 
Mot produce large variations in the predicted temperatures. 
Figures 20 and 21 show the actual temperature dis- 
uubutrons on the surface of the storage container and on 
Buessurface of the rocket motor respectively at maximum 
bulk temperature compared to a two dimensional TRUMP program. 
The TRUMP variation used for this comparison assumed no free 
Pemvection in the air gap and used tabular data to approxi- 


este the surface temperature of the storage container. 


C. GENERAL 
A comparison was made between surface temperatures on 

Per storage container that contained the rocket motor and 
the storage container that was empty. The low temperature 
was about 4°F colder in the empty container, whereas the 
Deren temperature was about 4°F higher on the container that 
contained the rocket motor. The empty container had a 
faster’ response time than the one containing the motor. 

"πο ditferences in heat capacities, eon erfeeis,.and 


Matura convection all contribute to the changes in tempera- 


ture noted. 


DD 








ca DE temperature Distribution at Surface of Storage 
Container at Maximum Bulk Temperature at approximately 
Ion 2 August 1972, 
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στε 213 Temperature Distribution at Surface of. the 
Rocket Motor at Maximum Bulk Temperature at approxi- 
mately 1500 on 2 August 1972, 
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P LIQUID CRYSTALS 

Meene apsulated cholesteric liquid crystals applied to 
Ehe surface of the storage container gave brilliant colors 
under the intense desert sun. These colors were much 
clearer and brighter than the same crystals viewed under 
laboratory lighting conditions. The liquid crystals photo- 
graphed well in both the color movies and the color slides. 
The movies showed by time lapse photography the rapidly 
changing surface temperature of the storage container. Two 
sample color prints made from the color slides are enclosed 
as Figures 22 and 23 to show the brilliance of the colors 
e reactor lity of obtaining data from color photos. 
The only photographic problem encountered was the intense 
Nono the Sunlight from the polyurethane film. This 
problem was partially overcome by taking the photographs 
from angles where the reflection was less intense. us br- 
tatively the liquid crystals were not adversely affected by 
the sun's rays after two weeks of desert exposure. No 
Beet sauantıitıtıve determination was attempted; however, 
rough approximations were made at the site. These approxi- 
mations were made by noting the color exhibited by a crystal 
at a certain time and then comparing the calibration of the 
crystal (Table 1) to the temperature recorded by the thermo- 
coup lettocated directly beneath the region of color change. 
The readings were within + 2°F, which was very encouraging, 
especially considering the approximations made while taking 


these measurements. Although photos were taken only during 
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the initial two weeks of the study, on site observations in- 
ica te that the crystals are still showing brilliant colors 
Beer 3-1/2 months. Preliminary RE indicates that the 
Eurethane film did protect the crystals from decomposing 
Ion the sun's rays and from being worn away by the wind 
blown sand. 

IS noted that the surface temperature of the storage 
container under the liquid cyrstals reached temperatures up 
momo F higher than a similar point not under the crystals. 
This 15°F difference was only evident when maximum tempera- 
Mies were obtained. During sunlight hours the temperature 
under the liquid crystals was always somewhat higher; how- 
Nat night both temperatures were equal. The difference 
in the container surface temperatures led to a difference of 
On. the surface of the rocket motor and 1°F at the center 
pane motor. It is believed that the difference in emissi- 
Sees of the gray and black surfaces resulted in the differ- 


ο -- in container surface temperature. 
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VI. CONCLUSIONS 


lou hecresudtscof this investigation, the following 
conclusions were dravn: 

NE Ithough-a sine wave is not a perfect fit for the 
EXE mental data at all points, it is useful in predicting 
IAN temperatures in the rocket motor, especially if only 
Aten h and low bulk temperatures are of concern. This 15 
especially true in the one dimensional case. In the two 
dimensional case, the problem of phase shift variations 
make the method of sinusoidal variation less desirable, 
EXUough still useful. 

Te simulation of the actual data by a table of 
temperatures gave the most accurate predictions of the 
experimental data. This method should be used whenever tabu- 
meda data are available; however, this will generally not 
Meche case for design work, in which case the sinusoidal 
approximations must be used. 

es lexibility ot both the analytical and computer 
simulations allow the changing of many parameters. The re- 
reine reffeecets of these changes on rocket motor temperature 
can be studied with the models, . | 

4. The convection assumption for this system resulted 
in only small changes in temperature and can be neglected 
When predicting design temperatures. Either the Liu or 
Beckmann correlation should be used to determine if convec- 


tion can be neglected in a particular system. 


62 





, Ι τ-ΕἘ otrtan empty storage container to obtain sur- 
face temperature data is a good approximation to using one 
fein a rocket motor inside. 

t is feasible to use liquid crystals for thermal 
mapping under desert conditions. Color photography with 
meaamdard equipment gives excellent results since brilliant 
SO Tors were observed. 

the liquid crystals appear to be stable for at least 
two weeks under the desert conditions when protected with a 
memsyurethane coating. 

de application of the liquid crystal system to the 
ace of the storage container resulted in large increases 
Ehe surface temperature of the container throughout the 
hottest part of the day. Care must be taken in applying 
and interpreting thermal readings from liquid crystals when 


exposed to radiant heating. 
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VIRT RECOMMENDATIONS 


Prom the results of this basic study, the following 
E mendations gor future work are offered: 

INTO refine the results of this project, a second 
rocket motor storage container system should be instrumented 
with the following changes: 

quid crystals shewld not be applied to the 
system used as the experimental model. As steel is a good 
thermal conductor, axial conduction on the surface of the 
Storage container may be significant. Heat flow from the 
area where the crystals are applied may lead to higher than 
Normal temperatures at other points on the surface of the 
container. 

b. The rocket motor should be weighed before and 
Beer the loading of the dry sand so that an accurate deter- 
mination of the density of the propellant simulant can be 
determined. 

c. Four additional thermocouples should be located 
ethe surface of the storage container to better enable 
the averaging of data. At present, the #1 thermocouple 
which was used as the average temperature reading of the 
top quarter of the surface of the container, in actuality is 
its hottest point; likewise the #10 thermocouple was used 
as the average temperature of the bottom quarter of the 
Surface, in actuality it's the coldest point. For averaging 


data, it is recommended that thermocouples be placed at 0130, 
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MO 0730, and 1030 and the quarters of the system be di- 
ed. at 0300, 0600, 0900, and 1200 to give a more realistic 
Dark temperature. Thermocouples at 1200 and 0600 will pro- 
vide the maximum and minimum temperature of the system. 

2. The TRUMP program should be rerun in both the one 
and two dimensional form, varying the mesh sizes to deter- 
mine the optimum number of nodes. 

5. A long term study of the effects of the desert en- 
vironment on liquid crystals should be done. The crystals 
Should be calibrated before Beine placed in the desert and 
Amen” brought to a laboratory for recalibration at specific 
intervals. 

4. Several modifications should be made to the TRUMP 
program to make it comparable to the version used at Law- 
rence Radiation Laboratory. The variable conductivity sec- 
won (BLOCK 2) and the PLOT subroutine (BLOCK 11) need to 
Memscorrected. The TIMEP subroutine which allows the setting 
of the problem time interval between data output should be 
added to this version of TRUMP. It would also be advanta- 
meous to increase the amount of tabular data that could be 
read in as boundary temperatures. 

5. From an academic standpoint, the effects of free 
convection in an air gap with varying boundary temperatures 


should be investigated. 
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APP DT xX A 


its came clon to Liquid Crystals 


Liquid crystals were first discovered in 1889 by Reinit- 
zer [Ref. 8] and the investigations of Lehmann which con- 
υεά tO 1915. Liquid crystals were considered to be 
oratory curiosities with no scientific or practical merit 
until the 1950's. They share some of the properties of both 
ias and crystals; for example, a typical liquid crystal 
substance scatters light in symmetrical patterns and re- 
flects different colors depending on the angle from which 
it is viewed. Studies in the last few years have helped to 
eny the unusual molecular structure of liquid crystals. 
Many applications arise from their ability to detect minute 
emeena ions in temperature, mechanical stress, electromag- 
fee radiation and chemical environment by changes in their 
eolor. 

Liquid crystals are divided into three classes; smectic, 
nematic, and cholesteric, depending on the degree of spatial 
arrangement of the molecules in the mass of the material and 
Mme type of the material {Ref. 9]. In this project only 
cholesteric liquid crystals were used and therefore only 
their properties will be mentioned. The molecular structure 
prscholesterie Liquid crystals is characteristic of the 
esters of cholesterol (Figure 24). The molecular layers 
are very thin with the long axis of the molecules parallel 


A ae ehe layers. The individual molecules are 
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ere 24: Molecular Structure of Cholesteric Ester 
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geure 25. Light Reflection from Liquid Crystals. 
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Uucucally flat, with a side chain of methyl groups (-CH,) 
projecting upward from the plane of each molecule. This 
EUUEI-curation causes the direction of the long axis of the 
molecules in each layer to be displaced slightly from the 
corresponding direction in adjacent layers, This displace- 
Ment, which averages about fifteen minutes of arc per layer, 
is cumulative through successive layers, and the overall 
displacement traces out a helical path. 

Mes molecutar Structure of cholesteric liquid crystals 
EM orise to many peculiar optical properties. If linearly 
polarized light is transmitted perpendicularly to the molec- 
mur c layers, the direction of the electric vector of the 
E11] be rotated to the left in a helical path. There- 
Hore, the plane of polarization will also be rotated to the 
tenes through an angle proportional to the thickness of the 
transmitting material. Liquid crystals are the most οσο μπι 
callv active substances known. Another strictly crystalline 
1151 property exhibited by cholesteric liquid crystals 
Mee ireu lar dichroism. When ordinary white light is inci- 
dent to a cholesteric material, the light is separated into 
Bro components, one with the electric vector rotating clock- 
wise, the other rotating counterclockwise. Depending on 
the material, one of these components is transmitted, and 
mme other is reflected. It is this property that gives the 
cholesteric phase its iridescent color when it is illumina- 
Hector jene. The particular combination of colors 
depends on the material, the temperature, and the angle of 


Encwimerdent light. 
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Ἱ πο ο οτ structure of a cholesteric substance is 
ver y delicately balanced and is easily upset. Any small 
disturbance that interferes with the weak forces between the 
molecules can produce marked changes in optical properties 
Esas reflection, transmission, birefringence, circular 
EM Drorsm; optical activity and color. The most striking 
«αἱ transformation that occurs in a cholesteric sub- 
Stance, in response to small changes in its environment, is 
Br ariation of color with temperature. The crystal lat- 
tice is disrupted by the thermal vibrations giving succes- 
sive transitions between the solid, the mesophase, and the 
meoeropic liquid with rising temperature. The change from 
Bmeomenree dimensional order of the crystal lattice to the 
disorder of the isotropic liquid occurs via one or more in- 
mediate states, each of which has a particular tempera- 
Mre range at which it is stable [Ref. 10]. 

À cholesteric liquid crystal system responds to changes 
in temperature by sequentially passing through the complete 
fetal spectrum (red through violet) in fractions or multi- 
degrees, depending on which cholesterol esters comprise the 
formulation. This color phenomenon is reversible and has 
been reported to function over a temperature range of -20°C 
Nc ery inportant point to note is that at a cer- 
tain temperature a given material or combination of materials 
will always exhibit the same color. Also, the rate of change 
from color to color as well as the exact temperature at 


voden the specific color changes occur are invariable. By 
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Ens cholesteric substances in various proportions, any 
desired temperature combination can be obtained. The thick- 
ness of the cholesteric film does not affect the predominant 
END" cnrz»th of the reflected light; the light becomes cir- 
muudarliy polarized [Ref. 11]. 

The colors scattered TM^ liquid crystals represent 
mus a fraction of the incident light (Figure 25). The re- 
maining portion of the incident light is transmitted by the 
liquid crystals. Therefore, an absorptive black background 
Dust be used to prevent reflection of the transmitted light, 
Bnereby enhancing the resolution of the scattered colors or 
wavelengths reflected by the liquid crystal system. 

Biezcholesterie Liquid crystal systems often present a 
number of problems due to the fact that they are viscous 
liquids. Some problems associated with the handling and 
the use of these materials are: 

sine tendency of the liquid crystal system to flow 
Pierre application can cause variations in applied film 
thickness. This may result in non-uniform thermal patterns. 

2. -Direct exposure of liquid crystals to adverse en- 
omental effects can cause variations in their sensiti- 
vity and deteriorate their color response in a few days. 

These problems can be partially overcome by using an 
encapsulated liquid crystal material system. The capsules 
are 20-30 microns in diameter and are a water-based slurry 
suitable for application by conventional coating techniques 


such as brushing or spraying. 
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Pmecapoutaved liquid crystals offer several advantages: 

mM They convert the liquid crystal system to a pseudo- 
EN. which provides for easier handling, application, and 
use. 

PEE hey provide longer shelf life by minimizing surface 
et amination and giving protection from ultraviolet light 
mef. 12]. 

hey exhibit relatively unlimited fatigue life. 

4. They reduce the ee dependence of the color 


observed, 


7k 





APPENDIX B 


Analytical Solution 


The method of complex temperature as presented by Arpaci 
Ber, 1] was used to find the steady periodic solution of a 
body experiencing a periodic sinusoidal disturbance. The 
general heat conduction equation in cylindrical coordinates 
was the basis for this derivation. It vas assumed that no 
heat sources existed in this problem, that the rocket motor 
storage container system was infinitely long, that there was 
RO heat conduction in the axial or circumferential direc- 
tions, and that the container surface temperature was spa- 

Eu gv uniform. Figure 26 gives a basic sketch of the system. 


The assumptions reduced the heat conduction equation to 


1 ο(νγ ---) 


Br 
r ΟΥ ae ot (1) 


Where r is the radial distance from the center of the rocket 
Motor, T is the temperature of the rocket motor at time t 
and position r, and @ is the thermal diffusivity, a property 
Πε conducting material. 


k 
uS (2) 


- 


ο το k is the thermal conductivity of the conducting 
Mea talas the density of the material, and c is the 
specific heat. All thermal properties were assumed to be 


constant over the temperature range of this problem. 
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Eure 26. Analytical Model of Experimental System. 
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The boundary conditions used in this derivation were 


and sz = - z (T - o) . at r = ro 


where τρ ΙΕ, ΠΠ; ΕΠΗ of the rocket motor. 
TI is the known storage container temperature which is 


assumed to vary as 


Den Tu” Τι) sin wt + T, 
where 
Tu = maximum bulk temperature of the storage container 
Pa = average bulk temperature of the storage container 
s : m ; 27 
w = frequency of the sinusoidal variation (ση mem 
t = time 


ü is the effective heat transfer coefficient across the 
air gap between the storage container and the rocket motor. 
It combines the heat transfer effects of radiation, convec- 
Bons and conduction into one coefficient. The radiation 
efficient was linearized by assuming constant temperatures 
(Ti, Το), representative of the average temperatures expec- 


ted in the problem, in the equation 


2 2 
heap ~ 71-2908, * T3, + T3“) 
where O is the Stefan-Boltzmann constant and on is the 
radiation exchange factor between surfaces 1 and 2, The 


convection coefficient is 


^ 


puce 
Doon 6 


where κ. is the effective conductivity of air as obtained 
from the Beckmann and Liu correlations [Ref. 5 and 7] and 


ο "το πιο of the air gap. In the analytical model, 
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mime effective conductivity was assumed to equal the conducti- 
ES thereby treating it as pure conduction and 


h = Baap * "CON 


Equation (1) was non-dimensionalized using the following 


relationships 


0 = pen (a non-dimensional temperature) 
D oe 
M 
E = τ (a non-dimensional distance) 


o 


to give 


TES τ 





mie 
i 

| 
e 
=| 


with boundary conditions 


go 0 at & = 0 


“and oe =-ß(8 - sin wt)at € = 1 


hr 





where 8 = is the Biot modulus (which compares the 

relative magnitudes of the effective heat transfer coeffi- 

cient across the air gap and the internal conduction 
resistances to heat transfer). 

An initial condition was not specified as the only 
Poncern ο ο Cae steady state, periodic behavior. Fol- 
Mowing Arpaci [Ref. 1], a complex temperature vas defined as 

P(r,t) = 0 (r,t) + 10(r,t) 


where Y(r,t) satisfied 





aces) Υ 2 
1 t BEEN 
PW i > 
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meen boundary conditions 


20.99 «199-0 at E = 0 


and oy = Ba). ο τος wt)+ i{-ß(8-sin wt) }at ἔξ 


This leads to O(r,t) which satisfied 


2 
, aces Po 08 


t DE Ue 





with boundary conditions 





a = 0 at € = 0 

do _ = 
and dE = -8(0 - sin wt) at £ « 1 
also Cote which satisfied 

do. 2 

1 o CES) το de" 

E ντ. τ 
with boundary conditions 

& 

4 = 0 at € = O 

"gs g 
and dE = ~8(0 -cos wt) at € = 1 


A solution of the form 


VRE) = HTI T(t) 


was assumed, where for large values of time T(t) was 
iwt 
assumed to equal e > therefore, 


iut 


Ur OD = olre, (4) 


Equation (4) was then substituted into equation (3) 


" dese) I ior “ὁ 
ri 


= 0 (5) 
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with boundary conditions 























ab = = () 
dt 0 | a € 
dd _ £ 
and m -B(o-1) at E = 1 
Equation (5) was expanded to give 
d? 9 iur - 
αφ ν ἆ se H 2 -Ὁ (6) 
de? E & α 
iur * 
Nov, let Z - S E 
and substitute into equation (6) 
2 
e x ox OS EE (7) 
zd 
dZ 
with boundary conditions 
d 
So = 0 at Z = 0 
and 
2 
d β ior, 
Si = = = > (9-1) at Z = » 
iur 
ο 
y a 
The general solution of equation (7) is 
d = C4 I (2) + C, K (Z) (8) 
as given in Ref. 13 with 
Ju co a 
MM TNCS (52) i 
IZ) = zZ) + + u = ee ee 
ο 2 (21)? (OPE 
! c ey ! 
and de 1 2n 
1 (72) ri 1 
E DN) ) E BEL es 


n=1 


Now using the boundary condition 


ol 
ce = 0 at Z = 0 


Jy 





mmaediftfferentiating equation (8) yields 
à d (I, (2)) d (K, (Z)) 
de TZ ee az 


where 
ἀ(1 (2)) 
Mz = 0 at Z = 0 
and 
d (K, (Z)) 
E- +0 at Z = 0 
therefore C, ZnO 
and d = σι 1. (2) (9) 


Now using the second boundary condition that 


p 2 
fiwr 
f 2 
CU 


O 
Ol 


and differentiating equation (9) gives 
d(1 (Z)) 
ἀφ = C 2 
dZ 1 az 


(1 (2)) 
Notane that BSc - 1,(2) and substituting into 


equation (10) 


q β pee | 
EUN e s τα /- 














1 
I  —————E—— 
ior ^ iur ^ ior * 
aß? οἱ α 
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and then substituting into equation (9) 














az) 
do = —— (11) 
te μα fu 
i α ) a BN ο ET \ α ) 
Now as 
J (imx) = i'I, (mx) fret 30 5. 135] 
2 2 
y Or 
I O jung) τν e O ¿3/2 ο 
Ο Q O Ol 
and 








2 2 
cor WT 
A ο το e E ( ο 3/2 
si o = - El α E e) 


Den 
Let a = W 2 amd substitute into equation (11) 








α 
Mar: at) 
φ = ————————————————————————— (12) 
ELDER 2 7372 τ 
J (i a) gi ' 8) 
3/2 ΠΟΙ 37 37 1 
Now i =e 4 = cos η + i sin πο σε (-1 + i) 
v 
Substituting this into equation (12) 
πο. AED 
pot (13) 
CARA 2 (121)5, (1/2 m 
O /28 
3/2 Deum 
AS J (aci ) = J, (ake 4) = BER (aé)+ iBEi (a) 


oud Neu 


J, (a£1?/?)- Ji (ate! 4) - BER, (at)* iBEi, (a£) 


Substituting these results into equation (13) yields 


BER (a£) + iBEi (al) 
a (14) 


΄ . a 
BER _(a)+ iBEi (a)* /3g (1-3) BER, (a) iBEi,(a)) 
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After rearrangement 


BER (a£) 4 iBEi (a£) 








m | 
[BER (a) ——BER, (a)-- ——B8Ei, (a) 1[BEi (a)t — BEi, (a)- —-BER, (a)] 
BB V2p Μπ y 28 
Letting | 
X = BER (a)+ —-BER_(a)+ --5-ΒΕ1. (α) 
R O /28 1 /28 T 
and 
X = 


BEi (a)+ EE) BER, (a) 
= /28 /28 


and substituting into equation (14) gives 


Piha jee i BE (αξ) 
ġ = yz AA 
Xp + iX, 


Rationalizing the denominator yields 


BER, (a£) + iBEi, (a£) 
αι ee δω. Q0») 


Xp sta X, 


Now 
E. (BER (a£) X, FBEi (a£)X, )Fi(BEi (a£)X,-BER (a£) X4) 
2 2 


Xp d Xi 


which after rearrangement gives 


2 LE 
yu (at) + BEi_ (ač) EH 


b = (νου 
2 2 
Rn + X, 
where 
S -1 BEi (a5)X, - BER, (a£) Xi 
= Stain 


BER (a£) X, + BEi_(a8)X, 
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Substituting into equation (4) gives 
ο oan RE) x 
s O ο i(wt +6 ) 
EN = πο e 
X dX. 
R n 


which also equals 


BER ^(a£)*BEi (ἀξ) 


Dam t) -'! 22 — [eos (ot*6^) i πι δν 
Xp T X: 


As this problem was modeled as a sine wave, the imaginary 
Bart of (r,t) was used. 
2 ο 
BER ^(a£)*BEi (aE) 


u en sin Con (17) 


Ap T X: 


Beh is the analytical solution of infinitely long concen- 
tric cylinders experiencing a periodic sinusoidal temperature 
Euruation on its outermost surface when heat conduction is 
assumed to be radial only. 

In summary 


"-— — 
T-T ΗΟ ΤΗΣ 





Α . κ 
toe ce Ate ο ο sin (wt+ 6 ) 
M A ` X LN 
R l 
where 
w ehet euperature of a point r in the rocket motor at 
time t 
1; ο ο ο υπο enp erature of the storage container 
Tu maximum bulk temperature of the storage container 
w = frequency of the sinusoidal variation(21/24 hours) 
t = time 
n ; . : 
E = Ia dimensionless distance from the center of the 
O rocket motor 
ΓΞ distance to the surface of the rocket motor 
ον απο from the center of the rocket motor 
az 
a = d = conduction parameter 
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α 5 i: - thermal diffusivity 
ο 


p = density 








e thermal conductivity 
c = specific heat 
BER = real Bessel Function 
fo — imaginary Bessel Function 
a Ξ a : 
Xp = BER, (a) + BER, (a) + ——— BEi, (a) 
/28 428 
: a : a 
X; = BEi (a) + — BEi, (a) - BER, (a) 
hr VIDIC /2 8 
B = p Biot modulus 
E -1 BEi (ab) Xp - BER (a£)X, 
of = tan 


BER (a6) X, Hi BEi (a£)X; 


The following computer programs were used to investigate 
a wide variety of parameters. The outputs are samples of 


some of these parameter studies, 
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APPENDIX C 


TRUMP Solution 


ΙΙ aseenputer program for solving transient and 
steady-state temperature distributions in multidimensional 
systems. This program was developed in 1965 at the Lawrence 
ED rou Laboratory by A. L. Edwards [Ref. 2] for their 
CDC/3600 computer. The program was adapted to the Naval 
Postgraduate School IBM/360 Model 67 computer system in 
EUMD; C. Exbayrum [Ref. 3] from ‘a version used by the 
Bees Goodrich Corporation. 

TRUMP is a multipurpose program able to solve a wide 
poruety of problems involving flow in various kinds of 
potential fields such as heat flow in a temperature field. 
TRUMP allows the solution of general nonlinear parabolic 
para] differential equations both in steady-state and 
transient problems. Complex geometric configurations with 
multidimensional flow may be solved using various coordinate 
Systems. Initial conditions may vary with spatial position. 
Material properties, boundary conditions, and other problem 
Pirsmeters may vary with spatial position, time, or the 
primary dependent variable. | 

Input data are fed to TRUMP in "Block" form through its 
I2 nput data blocks. A complete description of each of 
these blocks is given in Ref. 2. A model of the problem 
must be constructed and data from this model read into 


TRUMP through the data blocks. 
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Two models were used to simulate the rocket motor stor- 
age container system and several variations of each model 
mere investigated. 

The first model assumed one dimensional heat transfer 
(radial) with the assumptions that the System was infinitely 
long and that the container surface temperature was spatial- 
ly uniform. The system was modeled as two infinitely long 
EUUcentric cylinders separated by a 2.94 inch air gap. The 
wee Cylinder was constructed of 4130 steel and was filled 
with dry wind blown sand. The thermal properties of the 
Materials used in the experimental system are given in 
Table II with units most easily compared to the actual data 
obtained from the system at China Lake. 


FABLE τι 
Thermal Properties of Materials 


Material Density Specific Heat ΠΙΠεσπα] Conductivity 
Sand 0.054861bm/in” 0.195BTU/lbm’F 0.00026BTU/min-in-°F 
Steel 0.2807 lbm/in? 0.109BTU/lbm?^F 0.364  BTU/min-in-?F 
Air 079000436 Ἡ ο ο Tbm F 0. 0000225 a 


The model was subdivided into volume elements or nodes 
with the representative nodal points given in Figure 27. 
Although the representative nodal point may be located any- 
where in the node or on the surface of the node, in transient 
problems it is usually located so that the lines connecting 
the nodal points are perpendicularly bisected by the con- 
nected area. This gives maximum accuracy. Two boundary 
conditions were given to the surface node. The first was a 


sinusoidal disturbance which closely modeled the actual 
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iene 27]. Location of Nodes for One Dimensional TRUMP 
Model, 
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Eve experimental data obtained at China Lake as given in 
Appendix D. The ec wave exhibited a maximum temperature 
EN and an average temperature of 104?^F. Its period 
E. 4 hours (l440 minutes). The second boundary condition 
was the actual average surface temperature of the storage 
ο ἵσιποτ given at two hour intervals. Both these boundary 
BonWUrTtjons are approximations of the actual surface tempera- 
ture. Two hour intervals vere the minimum allovable for the 
tabular data as this version of TRUMP has a maximum table 
we of 12. 

Several assumptions were initially made. The thermo- 
Some data obtained from the experiment at China Lake gave 
the average temperature at a point on the storage container 
ΜΕ the actual outside surface temperature. As this 
Smmeainer wall was only 1/16 of an inch thick and made of a 
eood thermal conduetor, it was decided to model this data as 
a zero volume boundary node with a known temperature impressed 
A E as also assumed that heat transfer across the 
Em ap occurred only by radiation and conduction, neglecting 
the effects of free convection, 

It was estimated that the surface emissivities for the 
ο ο ταν and ele storage container were 0.9 [Ref, 6) 
basou on their haze gray surfaces. The radiation exchange 


factor for this geometry is given by 


PM ee 
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Sn le anput deck £61 the tabular approximation of the 
poundary condition is given at the end of this appendix. 
Eerat cycles of output data for the one dimensional model 
are also given., 

The second model assumed two dimensional heat transfer 
Goal and circumferential) with the assumption that the 
SS tem was infinitely long. The same physical model was 
assumed for the system except 48 nodes were used instead of 
IE Ihe representative nodal points are given in Figure 28. 
NI cur surface nodes (12, 24, 36, and 48) each had two 
different temperature approximations applied, a sinusoidal 
representation and a tabular input taken at two hour inter- 
vals. The four surface nodes were also modeled as zero 
volume boundary nodes. Each internal thermal connection 
between nodes is described in the input data by specifying 
HUE Vo node Identification numbers, two connector lengths, 
and two interface dimensional factors. An example of the 
thermal connections er node 4 is shown in Figure 28 and the 
πας data in BLOCK 5 of the two dimensional TRUMP program. 

Phe calculation of the radiosities in the two dimension- 
al case was accomplished by using a radiation-network and 
the E hod of crossed-strings. a 

The radiation shape factors for the two dimensional sys- 
tem were determined by the method of crossed-strings [Ref. 
''' τη: praphical construction for this method is given 


pure 29. 
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ντε «8: Location of Nodes for Two Dimensional TRUMP 
Model. 9 
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Figure 29: Graphi j 
μη ape Construction for Crossed Strings 
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I. is defined as the fraction of energy leaving sur- 

Mco vvhich directly reaches surface n. From the physical 
. . : : = == = sq». in = = 

dimensions of the model Ay A, Ας Qoa m A, Az As 

= 14,05 LL a, = 12,65 IL, same depth Let 

in Z in 

Ss equal the length of Αι. 

The crossed-string method lets each surface represent 
M@reneifective surface obtained by stretching a string tight- 
ly over the radiating face between the bounding edges, to 
produce a surface that cannot see any of itself. For ex- 
Aa surface 2' in Figure 29 stretched over surface 2. 


By definition F πας by reciprocity leads to 


2052 
Ay 
ea ok ea 
2 
and therefore since 
ΠΡ + Foot = 1 then Foo = 0.10 


lonc aTceulate the direct radiant heat exchanged between 
Surfaces l and 2, a minimum-length line was stretched con- 
to edge E of A, and a second minimum 


to edge F of A 


necting edge B of Ay 


Menta line from edge L of A These lines 


1 2 


are labeled Li υπ πι υ ο 559 and are equal to the width of 


the air gap, Li A ininum length lines were also 
Semetched from point B on Ay to F on A, and L on Ay to E on 
A). Ine length of these lines is Di and is made up of two 


parts; 51 the tangential distance from F to surface Ay 


πν,, 11 arc length from the point the tangent hits A1 


l 
nom E zs " 
to B. From geometry Di = Ys ro 6.62 
ο τς Ὁ 4.42" 
1 o a 


106 





therefore 


D c cp 11.04" 


1 1 1 
Zee 2L 
E 1 5 oz 
Now F1» =- x =. 0,86 
n 
meom reciprocity, Ay Pio = A, Foy 


Now Fi3 wS calculated from 
F1» + 2F,3 = ] 
therefore Fi4 = 0,07 


From symmetry δν ο Pis 


wma nen by reciprocity 


r 
- ο x 
Poy B Pr 0.047 
ih 
Now F was calculated by stretching minimum length lines 


23 


Mco E, from E to G, from F to G and from E to E; where 


the length of the line from F to E = from E to G = 5) 
from F to G = 2510 and from E to E = O 
2S,' - 2D 
2 1 
2 
Foc was που calculated from 
ν᾿... Ὁ Πε mn 
therefore Fos = 0.002 
As 


Fos was much smaller than the other ο it was not 
included in the radiation-network diagram in Figure 30 


Res 15]. 


1:07 








Eo" 








ο 


Figure 30: Radiation Network. 
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o calculate F ; 


ΠΕ the blackbody 


set to unity and all other blackbody 


potential of area l 15 


t 


potentials are set as 


zero. An energy balance was written at each node giving a 
set of six simultaneous equations as follows: 

Node 1 

E, Ay (1-3,) ‚= AQF15(J47J5)* 2AFy3 (4,733) 

Node 2 

EA, (0-3,) = A Fia (J359) 4 2A,F,, (J357J,)** 2A,F93 (37713) 

Node 3 

EJA¿ (0-32) = ALF, 5 (34-5, )+ AF, 4(J,-J5)+ AVF, 365,-J,)+ 


AF (I¿ IR AY) (3275) 


Node 4 

EA, (0-3,) = AJF)¿(I¿7I2)+ AF (3,7334 AJF),¿(J,735) 
Node 5 

EzA, (0-35) * 2A,F,4(J,-J,) 42A, F5  (J5-J,) FA F4 5 (15-4) 
Node 6 

EA, (0-J,) = 2AQF,,(Jg-J4)- AF, 5 (J5-J.) 


where > = radiosity of node n 








ı 
E^ - l-e a 
1 
E 
2 
o E 2 
2 
ΠΠ = radiation shape factors previously calculated. 


Now the values of A A mcd Fon were substituted into the 


energy balance equations which were then put into matrix 


foim 2s shown in Table III. 
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Betting Ere sF aan dando computer solution for matrix 


problems gave the radiosities as listed in Table IV. 


FABLE τν 
Radiosities at Nodes 


«9046 
70929 
.00495 
2000797 
‚000125 
«000080 


It i it ii ii 


αι C Q QûQ Q «€ 
ο Uu A U N E 


Iu o find the radiation exchange factors from node l to 
Nodes n, the radiation network shown in Figure 30 was re- 
duced to the equivalent network shown in Figure 31. Where 
the nodal equations are 


+ 


A, (1-0) = EA) (3,70) 


1-2 
where ROA 
τ 2 2 J = 71 
1-2 A, 2 
134, (1-0) -- E,A, (5,79) 
E, A 
2 2 - 
ID. a 33 7 «066 


These values of the radiation exchange factor are used 
in the two-dimensional program. A sample input deck for 
the sinusoidal boundary condition is included at the end of 
τοι... Several cycles of output data for the tabu- 


lar boundary condition are also given. 


rl 








igure 31: Equivalent Radiation Network 
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APPEND IOOCD 


Experimental Data 


mince Gata presented in this appendix were obtained from 
BH hermocouples on the rocket motor storage container sys- 
Bomelocated at China Lake, California. The thermocouple 
output was read out on a Honeywell Electronik 25, 24 channel 
E oder vhich had been calibrated at 50, 100 and 150°F. 
The data was taken on two consecutive, typical summer days 
Eust l and 2, 1972) at China Lake. Each thermocouple 
was read once every 24 minutes. The first set of data pre- 
sents the storage container tenperature at four locations 
ems three different ways of averaging this data. It also 
presents the ambient temperature and the approximate time of 
day. The second set of data presents the surface temperature 
of the rocket motor and three ways to average this data. It 
also presents the temperature at the center of the rocket 
motor and the approximate time of day. Figure 32 shows the 
Moeatdon of the thermocouples used to collect this tempera- 


fure data. 
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APPENDIX E 


Uncertainty Analysis 


oa certainty analysis was carried out on both the ana- 
Hytical solution and on a one dimensional TRUMP model of the 
rocket motor storage container system. In both models, the 
Smee tric heat capacity of the sand (9°), the conductivity 
of the sand (k), and the emissivity of the surfaces were each 
varied by ten percent to determine the sensitivity of the 
Sistem temperature response to each variation. Although 
Seem factors may also be varied, it was theorized that these 
amece had the greatest effect on the heat transfer of the 
system. These factors were also known with the least accu- 
racy; the maximum uncertainty of each was estimated to be 
plus or minus ten percent (odds 20 to 1). 

the analytical solution, varying the volumetric heat 
capacity changed parameter a, varying the emissivity changed 
Parameter B, and varying the conductivity changed both para- 
meters a and B. The effects on each parameter from each 


variation are given in Table V, 


TABLE V 
tinanoee in Parameters due to Changes in Thermal Properties 

Change in Property Change in Parameters 
Wollumetric Heat Capacity + 104 QN LP 

Volumetric Heat Capacity - 10% a mT 

Emissivity + 107 B + .49 

Emissivity - 10% B =- .37 

Conductivity + 102 a= ll, B= 22 
Conductivity - 10% αν ο ο ο Εις 2 
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Each factor was varied holding the other factors constant. 
The changes in temperature and time delay were computed from 
the difference between these new values and those previously 
med from the analytical solution. To obtain uncertainty 
mounds on the analytical curve, the second power equation 
[Ref. 16] was used, namely 


E 2 2 2 
W = yo, + Qu + p 


where 


Wp = resulting uncertainty in the calculated temperature 


due to uncertainties in temperature caused by 


ως estimated uncertainty in volumetric heat capacity 
me = estimated uncertainty in conductivity 
ως = estimated uncertainty in emissivity 


Eadentical calculation was carried out to calculate the 
uncertainty in time delay. The results of these calculations 
are shown in Figures 12 and 13 for the surface and center of 
the rocket motor respectively. The uncertainty in tempera- 
ture varied with time with a maximum variation of zo at 
Muclcenter of the motor and a maximum variation of +1.85°F at 
the surface of the rocket motor. The time delay varied by 
+31 minutes at the center of the motor and +1] minutes at the 
surface. The actual experimental data was also plotted on 
these Figures for comparison. 

The experimental data also had an uncertainty bound. 
Three primary factors made up this uncertainty bound; the 
accuracy of the thermocouple wire (+1.5°F), the readability 


ΓΕ πο recorder (+1°F), and the variation in temperature 
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caused by inaccuracy in the placement of the thermocouples 
(+1%F, estimated). The overall uncertainty in the experi- 


mental data was also calculated from the second power equa- 


Ron as 


= 2 2 2 = ο 
W ee een Zo 


οσο uncertainty bounds are also shown in Figures 12 and 13. 

procedure, similar to that used to find the uncertain- 
E oft the analytical solution, was used to analyze the 
Ne-UMEEng uncertainty in ae TRUMP numciical calceulation, 
M euts of these calculations are shown in Figures 14 
and 15. The uncertainty in temperature varied with time 
with a maximum variation of +2.95°F at the center of the 
rocket motor and a maximum variation of 41.95^F at the sur- 
XD i the motor. The time delay varied from +20 minutes 
Eu center of the motor to t9 minutes at the surface of 
ος motor. 

πα τος basis of the propagation of uncertainty analysis, 
it was determined that the solutions were most sensitive, in 
order of importance, to changes in the volumetric heat 


ποστς, emissivity, and the conductivity. 
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